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Employing some recent results in dynamics of systems with 
invariant subspaces we find evidence in both truncated and 
full axisymmetric mean-field dynamo models of a recently 
discovered type of intermittency, referred to as in-out inter- 
mittency. This is a generalised form of on-off intermittency 
that can occur in systems that are not skew products. As 
far as we are aware this is the first time detailed evidence 
has been produced for the occurrence of a particular form of 
intermittency for such deterministic PDE models and their 
truncations. The specific signatures of this form of intermit- 
tency make it possible in principle to look for such behaviour 
in solar and stellar observations. Also in view of its generality, 
this type of intermittency is likely to occur in other physical 
models with invariant subspaces. 



An important feature of the Sun and stars is their vari- 
ability on a wide range of time scales JjJ. In addition 
to the directly observed nearly periodic magnetic solar 
cycles, with an average period of 22 years, a particu- 
larly important feature of this variability is the presence 
of episodes, such as the so called Maunder Minimum of 
the 17 th century, during which solar activity (as deduced 
from sunspot numbers) virtually vanished [Q . The proxy 
data seems to indicate that this episode was not singular, 
but was preceded by numerous similar events occurring 
intermittently with an intermediate time scale of the or- 
der of 10 2 years. There is also some evidence suggesting 
similar variability in solar type stars || . 

Given the absence of naturally occurring mechanisms 
with such time scales in the Sun and stars, as well as evi- 
dence for the presence of non-linear phenomena in stellar 
and solar magnetic activity, one of the more plausible 
suggestions has been to associate this type of variability 
with some form of dynamical intermittency in the magne- 
tohydrodynamical dynamos operating in the solar/stellar 
interiors 

The complexity of the nonlinear partial differential 
equations (PDE) modelling these regimes has essentially 
led to three approaches, in turn employing: 
1. Low dimensional models which encode the main fea- 
tures of dynamo models and use, for example, normal 
form theory (see e.g. Q). These models, however, may 
not necessarily possess properties which we shall argue 
below are generic in axisymmetric dynamos. 



2. Low dimensional truncations fi|9|-|i"l| . 

3. Direct numerical integration of the full PDE models 
(e.g. @||). 

Models of type 1 have produced many useful insights 
into the nature of dynamo models. Models of type 2 
have been shown conclusively to be capable of producing 
a number of different types of intermittency, including 
crisis intermittency pH and type I intermittency p5| . 
Numerical studies of models of type 3 have produced a 
rich set of dynamical behaviours (see e.g . [^3 16 1), as well 
as intermittent types of behaviour |17|-|19||. The prob- 
lem has, however, been how to make precise the nature 
of these numerically produced dynamical modes of be- 
haviour. This is crucial for two reasons: firstly to be sure 
that these models are indeed capable of producing inter- 
mittent behaviour and secondly to use this information 
to characterise precisely their nature in order to compare 
their predicted signatures with observational data. 

The main difficulty has been the absence of an ap- 
propriate theoretical framework underlying such systems. 
Here, we shall use recent results in the transverse stability 
of attractors with invariant subspaces [p^|2C| l to demon- 
strate the presence in axisymmetric mean-field dynamo 
models of a recently discovered type of intermittency re- 
ferred to as in-out intermittency. 

Most studies of stellar dynamos have relied on mean- 
field theory, which is the approach we adopt here. The 
standard mean-field dynamo equation is given by 
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where B and u are the mean magnetic field and mean 
velocity respectively and the turbulent magnetic diffusiv- 
ity r) t and the coefficient a arise from the correlation of 
small scale turbulent velocities and magnetic fields 
In axisymmetric geometry, eq. (1) is solved by splitting 
the magnetic field into meridional and azimuthal compo- 
nents, B = Bp + B^, and expressing these components in 
terms of scalar field functions B p — V x A<j), B^, = B<fi. 

In the following we shall use a finite order truncation 
of the one dimensional version of equation ([!]) along with 
a time dependent form of a, obtained by using a spectral 
expansion, of the form: 
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where A n , B n and C„ derive from the spectral expansion 
of the magnetic field B and a respectively, J-, TL and Q 
are coefficients expressible in terms of m,n and I, N is 
the truncation order, D is the dynamo number and v is 
the Prandtl number (see |1| for details) . 

We first of all briefly discuss the main features neces- 
sary for the appearance of in-out intermittency [ p^|j20| 
and introduce the necessary vocabulary in the context of 
axisymmetric dynamo models. 

I. Invariant subspaces: In addition to simplifying the 
resulting equations and hence aiding the numerical inte- 
gration, the presence of symmetry forces dynamo models 
to possess invariant subspaces. For example, the trun- 
cated model (|J) with N = 4 is a 12-dimensional sys- 
tem of ordinary differential equations (ODE) with two 6- 
dimensional symmetric and antisymmetric invariant sub- 
spaces given by M s = {0, B x , 0, A 2 , 0, C 2 , 0, B 3 , 0, A 4 , 
0, C 4 } and M A = {Ax, 0, 0, 0, B 2 , C 2 , A 3 , 0, 0, 0, 
-B4, C4} respectively. Similarly the PDE model (1) pos- 
sesses invariant subspaces Ms = {B(9) = B(—8), A{9) = 
-A(-6)} and M A = {B{6) = -B(-0), A{6) = A(-6)}, 
where 8 is the latitude. As a result, a trajectory starting 
in either subspace remains in that subspace for all times. 

II. Non— skew product: If the dynamics can be written 
as an evolution in an invariant subspace forcing the trans- 
verse dynamics (skew product) we cannot have proper 
in-out intermittency but can have on-off intermittency. 
It is a generic property that we do not have a skew prod- 
uct. 

A related feature that seems to aid the appearance of 
in-out intermittency is the presence of non-normal pa- 
rameters that vary the system within the invariant sub- 
space as well as outside it. In the case of the truncated 
system (|2j), both the dynamo number D and the Prandtl 
number v are clearly non-normal parameters, as they en- 
ter the equations for A n and C n respectively and these 
variables in turn are a part of the invariant subspaces 
Ms and Ma- We note that generic parameters are non- 
normal. 

Recent studies of systems with these features [ fi5|]2"o| 
show the presence of a number of novel types of dynam- 
ical behaviour, including in-out intermittency. 

To characterise in-out intermittency, it is best to con- 
trast it with on-off intermittency j2^] as both types of 
intermittency occur in systems with invariant subspaces. 
On-off intermittency can occur as the result of an insta- 
bility of an attractor in an invariant subspace. It mani- 
fests itself as an attractor whose trajectories get arbitrar- 



ily close to an attractor for the system in the invariant 
subspace while makingoccasional large deviations away 
from it intermittently E^] . It can be modelled by a biased 
random walk of logarithmic distance from the invariant 



subspace 22 



We say an attractor A exhibits in-out intermittency to 
the invariant subspace M A (or Ms) , if the following are 
true @: 

1. The intersection Aq = AnM A (or A = AnM s ) is not 
necessarily a minimal attractor, i.e. there can be proper 
subsets of A that are attractors (for on-off intermittency 
Aq is assumed to be minimal) . This means that there can 
be different invariant sets in Aq associated with attrac- 
tion and repulsion transverse to Aq (hence the name in- 
out) and these growing and decaying phases come about 
through different mechanisms within Ma (or Ms). 
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FIG. 1. In-out intermittency in the ODE model (2) with 
N — 4 and parameter values D = 177.7 and v = 0.47. 
The energy and parity are given by E = E A + E s and 
P = (E s - E A )/E respectively, where E A and E s are the 
antisymmetric and symmetric parts of the magnetic field en- 
ergy with respect to the rotational equator ("antisymmetric" 
(P=-l) and "symmetric" (P=+l)). The top panel shows evo- 
lution of an initial condition in Ma and the other panels a 
nearby initial condition not in Ma- In these panels, we have 
taken a Poincare section at Ai=0, reducing the ODE to a 
map for clarity and comparison. 

2. The minimal attractors in the invariant subspace are 
not necessarily chaotic; they can be periodic orbits or 
equilibria. Furthermore, the trajectory remains close to 
this attractor during the moving away or "out" phases, 
with the important consequence that during these out 
phases the trajectory can shadow for example a periodic 
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orbit while drifting away at an exponential rate [£0| (see 
also H). 

3. The asymptotic scaling of the distribution of laminar 
phases in the in-out case can have two contributions: 



P n ~ an- 3 /V-^ + 7 e ( - 5n) =h+ J 3) 



(3) 



where a, /?, 7 and 8 are positive real constants depend- 
ing on the bias of the random walk modelling the "in" 
chain and the probability of leaking into the determin- 
istic "out" chain p0| ]. The term I\ is the one from bi- 
ased on-off intermittency. The extra term I2 can cause 
an identifiable shoulder to develop at large n which can 
help to statistically distinguish this type of intermittency 
from on-off intermittency. 

4. If the system has a skew-product structure, in-out 
intermittency reduces to on-off intermittency P0| . 

Because of the above considerations, we expect that in- 
out intermittency will be more generally visible in such 
dynamo models. Here we look at mean-field dynamo 
models and find in-out intermittency in such systems. 

Since the ODE models are more transparent, we first 
look at the truncated system (||) with N = 4. Fig. [j] 
shows an example of in-out intermittency in this system 
at parameter values D = 177.7 and v = 0.47. One can see 
clearly the "periodic" out phases (second panel), where 
the trajectory of the full system shadows the periodic 
orbit in the antisymmetric invariant subspace Ma (top 
panel) . Also we can clearly see the exponential growth of 
the amplitudes of the transverse variables through several 
orders of magnitude (lower panel). Furthermore, the cal- 
culated scaling of the distribution of the laminar phases, 
shown in Fig. is compatible with a curve obtained by 
fitting the parameters in the scaling law (H). 
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FIG. 2. Scaling of laminar phases for the model (2) with 
N = 4, D = 177.7 and v = 0.47, where the shoulder at large 
laminar phases (which is the influence of I2 and a character- 
istic of in-out) may be discerned. 

These signatures, namely the periodicity of the attrac- 
tor of the system restricted to the invariant submani- 



fold, the periodic locking and the exponential growth of 
the out phases and the compatibility with the scaling 
(^) clearly show the occurrence of in-out intermittency 
for the truncated dynamo systems. To show that this 
also happens in the full PDE axisymmetric mean-field 
dynamo models, we first of all recall that such systems 
possess the ingredients required for the presence of in-out 
intermittency, namely the existence of invariant subman- 
ifolds (symmetric or antisymmetric), non-skew product 
structure and non-normal parameters. 

Guided by recent results [jl7| , we integrated the mean- 
field dynamo equations using the code described in |l3| ] 
and implemented by JL9]. 
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FIG. 3. In-out intermittency in the axisymmetric PDE 
mean-field dynamo model (1). The parameters used were 
r = 0.4, C a = 1.942, Cn = -10 5 (C a Cn ~ D), f = 0.0, 
with the usual algebraic form of a = Qo/(l + B 2 ) (see [Jl7| for 
details of the parameters). To enhance visually the periodic 
locking we time sample the series in the two upper panels. 

Fig. ^| and Fig. ^ give examples of in-out intermittency 
in these PDE models with different algebraic forms of 
a. As can be seen this behaviour can occur with the 
invariant submanifold being either antisymmetric (Fig. 
||) or symmetric (Fig. ||). We have also obtained similar- 
results using a dynamic form of a of a similar type to 
that used in the truncated models. In addition to the 
presence of periodic behaviour in the system restricted 
to the invariant submanifold (top panel), these figures 
clearly show the presence of locking during the out phases 
(second panel) with an exponential growth of the energy 
of the transverse modes (bottom panel). This behaviour 
mirrors very closely the truncated model shown in Fig. 
|l| as well as that expected to occur from the theory [^0| . 
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Preliminary results show that the scaling for the PDE 
model is also compatible with (^) , but a clear verification 
requires a much longer integration. We shall return to 
a more detailed study of this issue in the future. The 
presence of the main features necessary for the occurrence 
of in-out intcrmittcncy in these models, as well as its 
presence in related truncated models, constitutes strong 
evidence for the occurrence of in-out intermittency in 
these PDE models. 
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FIG. 4. In-out intermittency in the axisymmetric PDE 
mean-field dynamo model (1). The parameters used were 
r = 0.4, C a = 1.5, Cn = -10 5 , / = 07, together with an 
algebraic form of a due to Kitchatinov J24| . The two upper 
panels are shown as in Fig. |^. 

We note that similar types of behaviour have also been 
seen in mean-field dynamo models with different topolo- 
gies |ll|^| thus lending further support to the suggestion 
that such behaviour is likely to occur in other settings. 

In summary, we have found strong evidence for the 
occurrence of in-out intermittency in both truncated and 
PDE axisymmetric mean-field dynamo models. Given 
the specific signatures of such intermittent regimes, this 
makes it in principle possible to test the presence of such 
behaviour in solar and stellar observations. 

We also note that since dynamical systems are generi- 
cally non-skew product with non-normal parameters, we 
expect this type of intermittency to be present in models 
of other physical systems with symmetry as well as their 
truncations. 
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